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ABSTRACT 

In the gravitational lens system Q2237+0305 the cruciform quasar image geometry is 
twisted by ten degrees by the lens effect of a bar in the lensing galaxy. This effect can 
be used to measure the mass of the bar. We construct a new lensing model for this 
system with a power-law elliptical bulge and a Ferrers bar. The observed ellipticity 
of the optical isophotes of the galaxy leads to a nearly isothermal elliptical profile 
for the bulge with a total quasar magnification of 16^|. We measure a bar mass of 
7.5 ± 1.5 x lO 8 hf 5 1 X0 in the region inside the quasar images. 

Key words: galaxies: fundamental parameters - galaxies: individual: 2237+0305 - 
galaxies: spiral - gravitational lensing. 



1 INTRODUCTION 

Gravitational lensing provides a unique way to weigh objects 
at cosmological distances without any assumption about the 
connection between light and dark matter. Since the discov- 
ery of the first gravitational lens (Walsh, Carswell & Wey- 
mann 1979) several gravitational lenses have been found and 
this method has been used many times to explore the mass 
distribution of galaxies. In this paper, we model the lens 
system Q2237+0305 in order to weigh the bar in the lensing 
galaxy. 

The quas ar Q2 237+0305 (z q = 1.695) was found by 
Huchra et al. ( [L985D at the centre of an SBb spiral galaxy 
(z — 0.0394) that is situated in the outskirts of the Pega- 
sus II cluster. The quasar was later resolved into four images 
that are situated arou nd the cor e of the galaxy within a ra- 
dius of one arcsecond (Yee 1985; Schneider et al. 1988). 

Two fundamentally different approaches have been used 
to model the lensing galaxy. One was to fit a parametric 
mass profile with several fr ee parameters to th e observed 
quasar image configuration (Kent & Falco 198S ) ; the other 
to use a model of the light distribution of the galaxy and 
t o fit for the mass-to- l ight ratio as th e single free parameter 
flSchneider et al. 1988; ; |Rix et al. 199^ ) . The former approach 



ences a light bundle from the quasar is small compared to the 
cosmological distances between observer, lens and source. 
The lens can therefore be treated as a mass sheet at the 
position of the galaxy. Since the galaxy disk of 2237+0305 
is inclined with respect to the sky, elliptical surface mass 
distributions must be used in the models of this system. 

Interestingly, the position angle, counted counterclock- 
wise from north, of the major axis of the elliptical lens mod- 
els found by Kent & Falco ( |l988| ) was about 67°. This is 
almost parallel to the axis through images C and D and just 
between the angle of the inclination axis of the galaxy (77°, 
Yee 1988) and the angle of the bar (39°, also Yee 1988). 
This situation is shown in figure This was also found by 
Kochanck (1991) and Wambsganss & Paczyhski (1994) who 



was nat urally much more procioo in the reproduction of the 



used simple circular mass distributions with an additional 
quadrupole perturbation. When fitted to the observed im- 
age geometry, the direction of the pe rturba tion turned out 
to be close to the one Kent & Falco ( 198§| ) found for their 
model major axis. M ore recent investigatio ns by Witt, Mao 
fc Sch echter ( |l99E| ), Kassiola & Kovner (|l995| ) and Witt 
fll996| ) obtained the same position angle for the perturba- 
tion or major model axis. 

On the other hand, the bar shows up promi nently in 
CCD images of the galax y. It has been no ted (Tyson & 



»h ao ^ ii ;™ ggc geometry due to the greatojamihei: of fr°° Gorenstein 1985| ; [Yee 1988; ; fritz et al. 1992Q that it might 



parameters. 

The length scale over which the lensing galaxy influ- 



contribute significantly to the lensing in the system - ini 
tially this was actually an aid to explain the lens ing effect 
when only the i mages A, B and C were known (Tyson & 



Gorenstein 1985) 
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The motivation for this paper is the idea that the ap- 
parent misalignment of the predicted model major axis and 
observed galaxy inclination axis as shown in figure [l] is due 
to the lensing influence of the bar. In section ^| we construct 
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Figure 1. Illustration of the image geo metry of Q2237+0305, 
motivated by figure 1 in Kent & Fal co (jl98q ). The images are 
labelled usin g the convention by Yee (1988) with positions from 
Crane et al. (L991). The relative areas of the circles correspond to 
the radio flux ratios from Falco et al. ( |l996| ). The position of the 
galaxy centre is indicated with a cross. The long arrows indicate 
the directions of the galaxy inclination axis and the bar, as well 
as the position angle of 67°. North is up, East to the left. 



and analyse a lensing model that includes the bar compo- 
nent and takes the observed position angle for the inclina- 
tion axis of the galaxy into account. Section ^ deals with 
the implications of this model for the bar. In section ^ we 
finally discuss our results. We use a cosmological model with 
H = 75h 7 5kms- 1 Mp C - 1 , Q, = 1 and A = 0. 



2 THEORETICAL MODEL 

The lensing model we construct has two components, each 
with several free parameters. In this section, we introduce 
the components and determine the values for these that pro- 
vide the best description of the observations. 



2.1 The bulge 



Yee (1985) identified three components in the inner part of 
the galaxy: bulge, disk and bar. In the lensing models for 
this system, bulge and disk have been represented by just 
one effective component since there is only limited observa- 
tional data from the quasar image and galaxy positions to 
constrain free parameters of the model. For simplicity, we 
call this composite compon ent 'b ulge'. 

Moreover, Kochanek (1991) and Wambsganss & Pa- 



czyhski (1994) found that the quasar image positions and the 
galaxy position in this system do not constrain the parame- 
ters even of simp le lens models. In particular, Wambsganss & 
Paczyhski (1994) showed that for a circular power-law mass 
distribution with an external shear there is a whole family 
of models that fit the observations; they discovered that for 



this family there is a linear relation between the magnitude 
of the external shear and the exponent of the mass profile 
for a vast range of exponents. The shear and the mass expo- 
nent are degenerate and one needs more information than 
only the positions of the quasar images and the galaxy to 
break this degeneracy. If one uses a two-component galaxy 
model the shear contributions from the two components will 
also be degenerate since the resulting shear is degenerate. 

One way to get around the shear degeneracy is to use 
an elliptical mass distribution instead of a circular profile. 
In this case, the ellipticity-parameter of the mass distribu- 
tion replaces the shear-parameter as a free parameter of the 
model. An analogous degeneracy in the ellipticity can then 
be broken by using the observed ellipticity of the isophotes 
of the galaxy. 

Generalising the approach by Wambsganss & Paczyhski 



(1994), we accordingly modelled the bulge with an elliptical 
power-law mass distribution with a major axis position an- 
gle of 77°. There is some disagreement in the literature on 
the value of this position angle (f or exa mple Rix, Schneider 
& Bahcall 1992). Fitte & Adam ( |l994| ) showed that this is 
because the position angle of elliptical isophotes is increas- 
ingly twisted towards the bar with increasing dista nce fr om 
the galaxy centre. The value we adopted from Yee ( 1988] ) is 
identical with the position angle determined from the galaxy 
continuum ma p within a radius of one arcsecond from the 
galaxy centre (fFitte & Adam 1994) where most of the lens- 



ing mass is situated. Let e be the elliptical parameter, so 
that the ratio b/a of minor and major axis of concentric 
elliptical shells is 



(1) 



It is useful to express surf ace mass densities in units of 
the critical lensing density ( Schneider et al. 1992j ) £ cr it = 



Z) s /47rG-DdI}ds, where D s , Dd and Dds are the angular 
size distances between observer and source, observer and 
deflector (lens), as well as deflector and source. The surface 
mass density k of a power-law elliptical mass distribution 
in units of £ C rit is given by 



" { "' Jh) = W 



(2) 



6i and 82 are the coordinates on the sky as measured in 
a coordinate system oriented with the observed major and 
minor axis of the bulge. 9 e is the elliptical radius 



+ 



(1 + e) 2 (1-e) 



(3) 



v is the power-law exponent of the elliptical mass distribu- 
tion and Eq is a constant. In the analysis of our results we 
use the ellipticity e = 1 — ^ = since this is the value 
that is usually used in the observations. In Appendix ^ de- 
flection potential and angles for this mass distribution are 
described. 



2.2 The bar 

The light distribution of bars has a well-defined elong ated 



shape, and is non-singul ar and centrally condensed ( Sell- 



wood fc Wilkinson 1993 ). It is not straightforward to de- 
termine the true form of bar mass distributions from this, 
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so that we have to assume a model. Very simple models 
with these pro perties of the bar light are the Ferrers profiles 
(Ferrers 1877 ). They were used in dynamical stu dies of bars 
( Freeman 1966a,b,c ; Martinet & de Zeeuw 1988; Sellwood & 
Wilkinf on since they can be treated analytically. To 

model the surface mass distribution of the bar of 2237+0305 
we used two-dimensional Ferrers profiles of the form 



k(0i 



M 1 - - ■£§■) 




if Zx a- Za. < i 
otherwise 



(4) 



In this equation, Q\ and 62 are the coordinates on the sky as 
measured in a coordinate system oriented with the observed 
major and minor axis of the bar. A is a real number, k c is 
the central surface density in units of the critical lensing 
density E cr it defined in section 2.1 and a, b are the semi- 



major respectively semi-minor axis of the bar. 

In the analysis we restricted ourselves to moderate ex- 
ponents A = 0.5, 1 and 2. The deflection potential and an- 
gles for integer values of A can be calculated analytically 
as described in Appendix 0. For A = 0.5 the profiles were 
constructed through the numerical superposition of many 
elliptical slices o f constant density (A = 0) and different size 
jSchramm 1993). 



2.3 The effect of shear 

The lensing influence of the bar can be understood by con- 
sidering a system with two shear tensors with shear direc- 
tions as shown in figure [l] for gal axy in clination axis and bar 
(see Schneider, Ehlers & Falco (1992) for the definition of 
the shear tensor) . The resulting shear tensor can be found by 
adding up the single tensors and the resulting shear direction 
is determined by the shear ratio of the two components. 

A source almost directly behind the core of the lensing 
galaxy appears lensed with four of the five images in a cross 
formation aligned with the axes parallel and perpendicu- 
lar to the resu lting shear direction, while the fifth is seen 
in the centre (Schneider et al. 1992, p252). In the case of 
2237+0305, the major axis position angle as found by the 
one-component lens models (67°) can be interpreted to be 
this resulting shear direction, which almost coincides with 
the axis through images C and D. The axis through images 
C and D has effectively been twisted away from the galaxy 
inclination axis by the bar. 

In figure ^ this twisting is illustrated by plotting the 
critical lines, caustics and image positions of two barred 
lenses with identical source positions. The caustics are the 
lines in the source plane that separate regions of different 
image multiplicity. The critical lines are the corresponding 
lines in th e lens plane where pai rs of images are created or 
destroyed ( {Schneider et al. 1992 ). 

In this figure, the first lens has a weak bar that barely 
changes the elliptical shape of the bulge's critical line or the 
corresponding diamond shape of the caustic. The other bar 
is significantly more massive; it warps the shape of these 
structures and shifts the image positions. 



2.4 Detailed modelling 

Our barred galaxy model has ten adjustable parameters. 
These are the positions of the galaxy and the source plus 
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Figure 2. Illustration of the effect of a strong bar. Assuming 
identical source positions, the diamond-shaped caustics, the cor- 
responding critical curves and the image positions are plotted for 
two different lenses. The smaller caustic and the almost elliptical 
critical curve belong to a model with a v = 1 power-law bulge 
and a A = 2 Ferrers bar that was fitted to the observed parame- 
ters of Q2237+0305. The positions of the images created by this 
model are labelled with unprimcd letters as in figure hi. In this 
model, the bulge mass inside the ring of images is about 20 times 
larger than the bar mass (see table hh. The caustic and the critical 
curve transform into the two elongated curves if the mass of the 
bar inside the ring of images is increased to half the mass of the 
bulge, while the bulge mass is kept fixed. The images are shifted 
by the more massive bar to the positions labelled with the primed 
letters. 



the constant Eq, the ellipticity e and the exponent v for the 
bulge as well as the bar mass normalisation k c , the semi- 
minor axis b and the exponent A for the bar. This number 
can be reduced by two if the observed ellipticity of the bulge 
and only fixed values for A are used. For the length of the 
semi-major axis of the bar we used the observed value of 
a ~ 9 arcsec (taken from the figures in Yee 1988 or Irwin 
et al. 1989). The lens effect is insensitive to the precise value 
of a because a is much larger than the radius of the ring 
of images (~ 1 arcsec). The length of the semi-minor axis 
must, however, remain a free parameter of the model since 
it is comparable to this radius, but not known well enough 
(b ~ 1 — 2 arcsec, see section ^|). 

There are ten observational constraints the system im- 
poses upon theoretical models. These are the coordinates of 
the four observed quasar images and the galaxy centre. The 
positions for the im ages and galaxy centre were taken from 
Crane et al. (1991). These positions have been determined 



from Hubble Space Telescope (HST) observations and have 
quoted measurement errors of 0'.'005. In general, the ratios 
of the fluxes of the different images of a gravitational lens 
also provide good constraints for a model. Unfortunately, in 
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t he ca se of Q2237+0305 t he ligh tcurves from Corrigan et al. 
(1991) or 0stensen et al. ( 1996| ) clearly show that all optical 
image fluxes are subject to flux variations due to microlens- 
ing of the quasar light from the stars in the lensing galaxy. 
In addition, the light from the quasar is non-uniformly dust 
reddened during the passage through the galaxy. The fluxes 
were, therefore, neglected in the modelling procedure. We 
will, however, compare the model predictions with the r e- 
cently measured radio flux densities by Falco et al. (1996). 



Let Ok , 6 g be the positions on the sky the model predicts 
for quasar images and the galaxy centre and O^o, go the 
observed positions with their positional uncertainties a^, a g . 
To find the best fit model, the expression 



4 

fe=i 



Ok — §ko) (8g — Sgo) 



+ ■ 



(5) 



(Wambsganss & Paczyriski 1994) was minimised through 



variation of the model parameters using a multidimen- 
sional minimisation routine (direction set or downhill sim- 
plex methods according to Press et al. 1992). 

In order to find an estimator for the separations Gk — dko 
between modelled and observed images for the first term on 
the right han d side of Equation (™, we used the method 
by Kochanek (1991); the separations between an optimally 
weighted source position and the positions in the source 
plane where the observed image positions are mapped to by 
a given lens model are propagated back into the lens plane. 



2.5 Analysis 

For a nearly circularly symmetric lens with a source almost 
in the orig in, it follows from Newton's theo rem in two di- 



mensions (Foltz et al. 1992; Schramm 1994) that the mass 



inside the circle of images is approximately given by the 
separation A8 of the images at opposite ends of the cross 



M = 



16 G D, 



DdDs (A8f 



(6) 



(see for example Narayan & Bartelmann 1996). The galaxy 
2237+0305 is situated relatively close to us at an angular 
size distance Z?d = 0.15 h^g 1 Gpc, so that D^ s /D s w 1. The 
separations of the quasar images are A9 « 1.8 arcsec, so that 
we get M w 1.5 x 10 10 h~^M Q . This value was also found 
in previous models for this system (Rix et al. 1992: 1.44 ± 
0.03 x 10 10 hf 5 1 A4 Q , Wambsganss & Paczyhski 1994: 1.48 ± 
0.01 x 10 10 Ii^Mq). The other value theoretical models for 
2237+0305 agreed on was the resulting shear direction of w 
67°. 

In our barred lens model, the bulge acts as the main 
lensing mass and the bar as a perturbation; for a given el- 
lipticity e or exponent v, the bulge parameter Eo and hence 
the bulge mass do not change very much for different A-bar 
models. In fact, experiments with different bar masses as in 
figure ^ showed that for similar masses of bulge and bar in- 
side the quasar images the cruciform image symmetry gets 
skewed. Models with a strong bar thus cannot reproduce a 
symmetric image geometry as in Q2237+0305. In order to 
explore the effect of A on the models, we used fixed values 
A = 0.5, 1 and 2. 




0.75 



0.5 e 



0.25 



Figure 3. Contour plot of confidence regions in the parameter 
space of bulge exponent v and ellipticity e. The indicated contours 
contain 68.3%, 95.4% and 99.7% of normally distributed models 
around the minimum of x 2 • 



Besides Eq and A, the parameter space of e, u, k c and b 
had to be examined. We first scanned the parameter space of 
e and v while leaving n c and b as free parameters. In figure ^ 
a contour plot of the confidence regions that contain 68.3%, 
95.4% and 99.7% of normally distri buted models around the 
minimum of x 2 ( Press et al. 1992 ) in the parameter space 
of v and e is shown. Only the plot for A = 2 is presented; 
the cases with A = 0.5 and A = 1 are very similar. The 
parameter space was scanned with a stepsize of 0.02 for v 
and 0.01 for e. For every point the best parameters were 
determined through minimisation and a x 2 - vam e was com- 
puted. The best fit models lie in a long valley that extends 
up to v ~ 1.25. The unclear structure at the lower end of 
the valley for v < 0.5 and e < 0.1 is due to numerical ef- 
fects. Detailed investigation shows that the valley continues 
towards smaller v, becoming shallower. The models in this 
region, with low e and v, are similar to circular disks with 
constant surface mass density and are not examined here 
because they do not represent realistic galaxy models. 

The external-shear models by Wambsganss & Paczyhski 



(1994) yielded no constraints on their circular mass distri- 
butions for the whole range of parameters from v = 0.07 to 
v — 2. The smaller allowed range for v in figure ^ shows 
that the ellipticity does not allow the same kind of freedom 
as the shear. For high ellipticities of the power-law bulges, 
the quasar image geometry cannot be reproduced anymore. 

In figure ^ the total mass, the bulge mass, and the bar 
mass inside a circle of 0.9 arcsec as well as the semi-minor 
axis b and the total magnification fi along the valley of best 
fits around the v = 1 profile are shown for 0.75 < v < 1.3 
and the three values of A. The plots for different A only differ 
noticeably in their semi-minor axis predictions. Beginning 
with v ~ 1.2 the minimisation produces numerical noise at 
the upper end of the valley from figure |^ since the fits get 
worse. 
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Figure 4. Model parameters along the valley of best fits. Plotted 
against v are the total mass and the masses of bulge and bar 
inside a circle of 0.9 arcsec in 10 10 bfg 1 Mq as well as the semi- 
minor bar axis b in arcsec and the total magnification fi. Different 
line styles have been used for different values of A as indicated in 
the b-v panel. 



It can be seen that the total mass of the model inside 
a circle of 0.9 arcsec is almost constant. The constituent 
masses of the bar and the bulge exhibit a dependence on v. It 
can be interpreted that the bar mass increases with v in or- 
der to counter the increased contribution of the bulge to the 
resulting shear due to the increase of e, and hence the bulge 
mass decreases in order to conserve the mass inside 0.9 arc- 
sec given by equation ^ . The magnification drops strongly 
with increasin g v, w hich was also observed by Wambsganss 
& Paczyriski ( 1994| ). 

As motivated in the introduction, we chose the model 
from the family of best models in figure ^| that exhibits the 
observed ellipticity o f the bulge. The ellipticity has been 
measured by Racine (1991) to be e = 0.31 ± 0.02, in agree- 



ment w ith th e results from the continuum map by Fitte & 
Adam (1994). This range of ellipticities is encompassed by 



the range v — 1.05 ±0.1 of bulge exponents. We use this re- 
gion of allowed bulge models to determine the uncertainties 
of the predictions of models with different bars. 

From all possible values for v and the three bar models 
we obtain an estimate of the total mass inside a circle of 0.9 
arcsec of 1.49 ± 0.01 x lO lo h~ 1 A4 . This is consistent with 
the estimate and the results from the literature given at the 
beginn ing of this section. In a review, Narayan & Bartel- 
mann ( 1996 ) mention only two cases of gravitational lenses 
in which it has been possible to constrain the radial distri- 
bution of the lens galaxy. Both models are based on singu- 
lar elliptical or nearly elliptical profiles. For MG 1654+134, 




-1.5 b 



Figure 5. Contour plot of confidence regions in the parameter 
space of bar normalisation k c and semi-minor axis b (in arcsec) for 
a bulge exponent v = 1.05, ellipticity e = 0.31 and bar exponent 
A = 2. The indicated contours contain 68.3%, 95.4% and 99.7% 
of normally distributed models around the minimum of \ 2 . 



Kochanek fll995| ) obtained v = 1 .0 ± 0.1, and for QSO 
0957+561, Grogin & Narayan ( [l996| ) obtained v = 1.1 + 0.1. 
The predictions of our model for the bulge are in good agree- 
ment with these values. 

In table ^, the model parameters for the different values 
of A and the uncertainties due to the uncertainty in v are 
shown. The x 2 -values from equation (^) have been divided 
by the number of degrees of freedom, giving a measure of the 
quality of the fit. The number of degrees of freedom is the 
number of constraints minus the number of free parameters; 
here we have three degrees of freedom since e, v and A are 
fixed. 

The value for the total magnification predicted from the 
v = 1.05 mod el is a bout half the value found by Wambsganss 
& Paczyriski ( 1994 ) for their v — 1 circular power law profile 
with an external shear. This illustrates that the use of an 
elliptical mass distribution drastically changes the predic- 
tions of the model; a similar discrepancy is apparent for the 
time-delays. Note that /i to tai » f = 16.7 with e = 0.18 (cor- 
responds to e = 0.31, see section 2.1) as derived for a v = 1 
power-law lens with smal l e and a source in the origin by 
Kassiola & Kovner (1993). In co ntras t to the non-singular 
mass models by Kent & Falco (1985) our lens models do 
not produce a fifth image in the centre due to the central 
singularity of the mass distribution. The predicted source 
positions are very similar to the ones by Kent & Falco. 

The remainingparameter space of the bar models is il- 
lustrated in figure g| Similar to figure |^, the contours of the 
confidence regions of normally distributed models around 
the minimum of x 2 in the parameter space of k c and b are 
shown for the model with v — 1.05, e = 0.31 and A = 2. 
The parameter space was scanned with a stepsize of 0.003 
in k c and 0.025 arcsec in b. The bar parameters are well con- 
strained at the bottom of a steep boomerang-shaped valley. 
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Table 1. Model Parameters for Q2237+0305 for three values of the bar exponent A and the bulge exponent 
v. The columns for A = 0.5, 1.0 and 2.0 are indicated. Each table entry contains the value for v = 1.05 
and the differences to the corresponding values for v = 1.15 (upper index) and 0.95 (lower index). Only 
one value is given for a line if the differences between the bar models are below the rounding precision, x 2 
is divided by three, the number of degrees of freedom. Eq, b and the source position (/3i,02) are given in 
arcsec. Mtotal is the total magnification, /xjj are the relative magnification ratios between the images. A^j are 
the relative time— delays between the images in h^g 1 hours. iWbulge(< 0.9") and M Dar (< 0.9") are the masses 
of bulge and bar inside a circle of 0.9 arcsec and M Dar total i s t ne total mass of the bar. Masses are given 
in 10 10 h^ 1 A^Q. The semi-major bar axis was a ssumed as a = 9 arcsec. The VLA flux ratios are the 3.6cm 
radio flux ratios measured by Falco et al. (1996). 
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4.9 



+0.4 
0.3 
+2.7 
4.4 
+ 1.0 
0.8 



A/ b uigc(< 0.9") 
M bar (< 0.9") 

-^bar, total 



1 42~ u ' uz 
1, ^ Z +0.01 

07+ - 01 
u - u '-0.01 

47+0-07 

u '^'-0.06 



1.42 
0.07 
0.43 



0.02 
+ 0.01 
+0.01 
0.01 
+0.06 
0.06 



1 41 -"-"J- 

o.08+°;£ 



Table 2. Local lensing parameters at the positions of the quasar 
images. The values have been determined with a A = 1 bar, but 
they are identical for A = 0.5 and 2 within ±0.01. Each table 
entry contains the value for v = 1.05 and the differences to the 
corresponding values for u = 1.15 (upper index) and 0.95 (lower 
index) . 



Image 



7 



A 
B 
C 
D 



0.36 
0.36 
0.69 
0.59 



+0.07 
0.06 
+0.07 
0.01 
+0.01 
0.03 
+0.02 



40+ 03 
u -* u -0.03 
,+0.03 
-0.03 
,+0.09 
-0.09 
,+0.06 
-0.06 



0.42Z 

0.7lj 
0.61"' 



The surface mass density k in units of the critical den- 
sity and the local shear 7 at the positions of the images are 
given in table ^. The values for v = 0.95 and 1.05 are sim- 
ilar to what one gets for a circular v = 1 power-law mass 
distribution (singular isother mal sphere) with a mixture of 
inte rnal and external shear ( Kochanek J.991 ; Witt fc Mao 
1994 |). ]| t can be taken from tables |l| and ^ that a change of A 
does not cause a measurable change of observable quantities 
except for the semi-minor axis of the bar; the uncertainty 



of the model parameters is dominated by the uncertainty of 
the bulge ellipticity e/exponent v. 



3 PROPERTIES OF THE BAR 

The width of the bar has not been meas ured previously. In 
figure 1 of their paper, Irwin et al. (1989) present a contour 
plot of the galaxy where the bar has been separated from 
the disk. In this plot, the bar appears about 18 arcsec long, 
but only « 2 — 4 arcsec wide. It is just the less well-known 
minor axis that enters the lensing model since the quasar 
images are situated in the centre of the galaxy. In their image 
analysis Irwin et al. subtracted structure from their image 
that is smooth on scales of ~ 5 — 10 arcsec. This procedure 
removed the galaxy disk very efficiently, but, being long as 
well as thin, the bar shown in their figure 1 could also be 
affected by this procedure. 

A different approach to obtain the li ght distribution of 
the bar was pursued by Schmidt (1996). Using a Hubble 



Space Telescope (HST) I-band image taken with t he Plane- 
tary Camera prior to the first servicing mission ( Westphal 
1992), the galaxy was decomposed into bulge, disk and bar 
with analytical profiles that have been convolved with the 
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point-spread function of the telescope. Exponential profiles 
(Andredakis & Sanders 1994) were fitted to bulge and disk. 



After these components were removed from the image, the 
bar and the spiral arms remained. The bar light distribution 
could be fit with a A = 2 Ferrers profile with a — 9.5 ± 1.0 
arcsec and b = 1.0±0.3 arcsec. This light model can be com- 
bined with the mass model for A = 2 from table [j] due to the 
similar value for b; the I-band mass-to-light ratios of these 
model components inside a circle of 0.9 arcsec are given by 
M/Li w 4.8 1175 for bulge plus disk and M/Li ^ 5.0 1175 for 
the bar. The bar mass detected with gravitational lensing 
and the bar light in the I-band in this model both consti- 
tute a 5% fraction of the total mass respectively light inside 
0.9 arcsec. 

This result is, however, dependent on the bu lge model. 
For a de Vaucouleurs-bulge (de Vaucouleurs 1948), a A = 0.5 
Ferrers profile fitted the bar light much better with a similar 
value for a, but a much larger value b — 3.1±0.9 arcsec which 
cannot be combined with the lensing model from table [l] 
since b is very different. The unrefurbished HST point spread 
function inhibited a clear distinction between the quality of 
fit of these different bulge models, so that the question about 
the minor bar axis is not decided yet. 

A value of b > 1 arcsec could in connection with the val- 
ues for b in table ^ be taken as evidence in favour of steeper 
bar models, A > 2. There is observational evidence from the 
light d istributions of real bars for more boxy mass distribu- 
tions ( freeman 199(: ), so that it has to be ascertained that 
this result is not a relic of the approximation of the bar with 
an elliptical shape. At the positions of the quasar images, 
however, an ellipse is a good approximation of a box due to 
the large semi-major axis of the bar, a « 9 arcsec. A lensing 
model with a steeper Ferrers profile with A — 10 leads to 
a larger bar with b = 2.0. This indicates that a measure- 
ment of the extent of the bar could in principle constrain 
the mass distribution of the bar. If one wants to use this 
method, high resolution images with small seeing have to 
be used since an accuracy of the bar width on the order of 
tenths of an arcsecond is needed. 



4 DISCUSSION 

In this paper we have presented a barred galaxy model for 
the gravitational lens 2237+0305. We used a power-law el- 
liptical mass distribution for the bulge and chose the model 
for which the observed ellipticity is predicted. It turned out 
that this model has an exponent close to unity, which is 
compatible with other determi nations of lens mass profiles 
through gravi tational lensing (Kochanck 1995; Grogin & 
Narayan 1996). The bar represents a small perturbation of 



the deflection field of the bulge of the galaxy, amounting to 
7.5 ± 1.5 x 10 s ib^MQ or about 5% of the bulg e mass m 
the critical region inside the quasar images. 

The relative magnifications our model predicts for the 
quasar images can be compare d with the 3.6cm radio flux 
ratios published by Falco et al. ( 1996[ ) . Their results are also 
given in table |l| Falco et al. argue that it is unlikely that 
the radio flux densities are variable from microlensing due 
to the larger size of the radio emitting region as compared to 
the optical continuum emitting region, although they cannot 
completely rule out microlensing as an important effect in 



the radio. If microlensing is not important, the model mag- 
nification ratios should be identical to these measured flux 
ratios. It can be seen that only the ratio /xda between im- 
ages D and A is not compatible with their results although 
no effort was made to fit the flux ratios. 

In order to find out more about the discrepancy of 
the ratio ^da between observation and model one has to 
make the relatively large error bars from Falco et al. smaller 
through longer radio observation of the object. Unfortu- 
nately, Q2237+03 05 has a radio flux density of only pb 1 mjy 
(Falco et al. 1996), so that radio observations of this object 



are very time-consuming; Falco et al. observed for 11 hours 
of which only five could be used eventually due to weather 
conditions. 

To get additional, independent arbiters for the model, 
it would be very helpful to measure the time-delays in this 
system. Since the time-delays are of the order of several 
hours, this has to be done in a wavelength domain where 
the necessary intra-day variability is likely to occur for a 
radio-quiet quasar, for example in t he x- ray regime as pro- 
posed by Wambsganss & Paczyhski (1994). Also, monitoring 
in the radio would show if the quasar image flux densities 
vary at these frequencies. Unless we learn more about the 
radio flux densities and the time delays, it is not possible 
to decide whether or not it is microlensing that causes the 
low magnification of image D. In the optical, image D has 
always been the faintest quasar image. In fact, in the first 
reso lved image of the quasar, image D was not visible at 
all ( Tyson fc Gorenstein 1985 ). There is also spectroscopic 
evidence from optical data that image D is undergoing de- 
magnification (Lewis et al. 1996). 

If image D is in fact microlensed in the radio, the con- 
sequences are interesting. The scale size of the radio region 
could be less than the characteristic scale of the caustic net- 
work. Alternatively, the radio source could have an asym- 
metric structure like a jet that would have differing mi- 
crolensing properties for different paths of the microlenses 
across the source. 

Yet another way to significantly change the radio flux 
density of image D would be a globular cluster or black hole 
(Lacey & Ostriker 1985) with a mass of about 10 6 7V(q in 
the halo of the lensing galaxy that is situated close to image 
D. An object of this mass would magnify or demagnify the 
radio image of the quasar, depending on its location with 
respect to the direction of the local shear. This effect, the 
perturbation of lens models by 10 s Mp^ o bjects, has recently 
been treated by Mao & Schneider (1997). With this, we can 



estimate that a surface mass density of globular clusters or 
black holes of approximately 0.04 E cr it (£ cr it is defined in 
section [2.l| ) or 470 1i75.Mq/pc 2 is needed to observe a de- 
magnification of image D by 40% or more with a probability 
of 20%. Higher surface mass densities would make it more 
likely. This is much more than the globular cluster surface 
mass density of about 1 Mc^ /x>c 2 seen in our Galax y within 
5 kpc of the Galactic centre (Mao & Schneider 1997). It thus 
seems unlikely that the demagnification is due to a globular 
cluster. 

The question of the existence of such a massive object 
near image D could be solved with a metho d that was pro- 
posed by Wambsganss & Paczyhski ( 1992 ). They showed 
that these objects would bend or even create holes in the 
radio maps of milliarcsecond jets of gravitationally lensed 
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quasars. If we could observe extended structure of the im- 
ages of the radio-weak Q2237+0305, features due to glob- 
ular cluster or black hole lensing could be easily identified 
because they would not be seen in the other gravitationally 
lensed images of the quasar. 

It has recently been suggested (Keeton, Kochanek & 
Seljak 1996; Witt & Mao 1997) that a strong lensing pertur- 
bation is necessary for a number of lenses in addition to an 
elliptical mass distribution in order to model the observed 
image geometry. In our model for Q2237+0305, the addi- 
tional perturbation from the bar is small, but the system 
is, nevertheless, unique in that we can see the perturbing 
agent. For the lens systems mentioned in these papers, the 
better fit for the models was obtained with an additional 
external shear. We saw, whe n we compared our results with 
Wambsganss & Paczyhski (1994), that the predictions for 
magnifications or time-delays from shear models differ by 
up to a factor of two from the predictions from elliptical 
mass deflectors. If one wants to get lensing models with reli- 
able predictions for magnifications or time-delays, the nec- 
essary perturbations will ultimately have to be generated 
i n the models by m ass components like dark matter haloes 
(Keeton et al. 1996) or bars. 
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APPENDIX A: POWER-LAW ELLIPTICAL 
MASS DISTRIBUTIONS 

The lensing properties of the elliptical power-law mass dis- 
tributions from Equatio n ([^) w ith v = 1 have been described 
by Kassiola fc Kov ner (1993) and Kormann, Schneider, & 
Bartelmann (1994). Kormann et al. determined these by 
solving the Poisson equation for the deflection potential ij) 
(see Schneider et al. 1992 for the definition of the deflection 
potential) in polar coordinates 



'2 and ip 



ld_ 

9 86 



M\ 1 d 2 ip „ ,„ . 



06 J 



(Al) 



Their approach can be generalised for arbitrary real numbers 
v with < v < 2. Define n = 2 — v, s n = sinn^, c n = cosn|- 
and 



Using the integrals 



(A2) 



A/i 



A/2 = 



cos nip 



la 1 luj . / 



A"(<pO 



sm nip 
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A/a = M + 



cos ntp 



7T°V, 



a* <y ) 

the deflection potential for < tp < ? is given by 



— A/i sin nip + 
n 



i(AA 2 + ^ 



A/3 



cos rw^ 



(A3) 



(A4) 



The potential for the other quadrants can be taken from 
this result because of the elliptical symmetry of the mass 
distribution. The deflection angles a. = Vi j) for this quad- 
rant can be calculated from Equation (A4). The cartesian 
components are 



Mi ( 



sm nip cos ip — cos nip sin ip 



+ ^A/2 H — — A/3^ ^ cos n<^) cos (^9 + sin sin ip^j 
02(9, ip) = Eq A/"i ^ sin ny> sin </p + cos nip cos 93^ 



+ (AA 2 + ^ 3 )( 



cos nip sm 95 — sm nip cos 95 



.(A5) 



These expressions can be evaluated numerically. For v — 
n — 1 the integrals are analytically solvable and the for- 
mulas become identical to t he res ults by Kassiola & Kovner 
(1993) and Kormann et al. (L994). After this wor k was com- 
pleted, we discovered that Grogin & Narayan (1996) also 
used power-law elliptical mass profiles in their model of the 
lens of the double quasar 0957+561. In their paper they 
present the deflection angles of power-law elliptical mass 
distributions in a complex- valued lensing formalism in terms 
of the complex hypergeometric function. 



APPENDIX B: FERRERS PROFILES 

The deflection potential ipo and the deflection angles c?o of 
an elliptical slice of constant surface density, a Ferrers profile 
from Equation (W) with A = 0, are given by 



ipo (01 



-^abK c (Q m - OiQio - OtQoi) (Bl) 



and the cartesian components 

«oi (61, 02) = abK c 6iQ w 
ao2 (61,02) = abK c 2 Qoi, 
where 



(B2) 



Qoo 
Q01 



= 2 In 



(Va 2 + P+ V b ' 2 +p) 



-b 2 



= 2/A(p)-. 



a 2 + p 
b 2 +p 



p is the positive solution of 

ei el 
+ .„ = 1 



(B3) 



(B4) 



a 2 + p ' b 2 + p 
outside the bar, and p = inside the bar. A (p) is given by 

A (p) = vV + /5 ) (b 2 + p). (B5) 



ao has first been derived by Schramm (1990) by using the 
known force field for a homogenous ellipsoid and letting the 



largest axis go to infinity in order to get the corresponding 
2-dimensional equations. To derive the potential ipo in terms 
of real- numbered coordinates 61,62, we app lied Schramm's 
method to the results by Pfenniger (1984) for the poten- 
tial and force fields of 3-dimensional Ferrers ellipsoids. This 
enables us to calculate deflection potential and deflection 
angles for a whole family of Ferrers profiles. The deflection 
potential for Ferrers surface-density profiles with integer ex- 
ponents A is given by 



ipx {61 



abK c 



2(A + 1) 



du 



A( 



a 2 + u 



b 2 



A+l 



(B6) 



The correspondi ng 3- dimensional expression was first de- 
rived by Ferrers (1877). For integers j, k, j 7^ or k 7^ 0, this 
integral can be split up into a sum containing the coefficients 



ijk = 



du 



A(w) (a 2 + u) J {b 2 + u) k 



(B7) 



For Qoo use Equation (B3); a remaining contributio n from 



the infinite axis had to be subtracted here. Pfenniger ( 1984 ) 
solved his corresponding integrals using recurrence relations. 
Translated into two dimensions, the Qjk obey the relation 



Qjk = (Qj-l,k — Qj, k- 

as well as for n > 



Qon 



2n-l 
1 



2n 



A(p) (a 2 + P r 
2 



A( P )(b 2 + P y 



(B8) 



(B9) 



(B10) 



With these equations, the deflection potential ip\ as well 
as the deflection angles Qx = Vip\ can be calculated from 



Equation (B6). The coefficients Qij can be treated as con- 



stant for the derivation with respect to 61 and 6 2 since the 
definition of p in Equation (B4) implies that = 0. 
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